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MEASUREMENTS  SHOWING  THE  FEASIBILITY  FOR  RADAR  DETECTION 
OF  HAZARDOUS  WIND  SHEAR  AT  AIRPORTS 

R.  B.  Chadwick,  K.  P.  Moran, 

G.  E.  Morrison,  and  W.  C.  Campbell 

NOAA/ERL/Wave  Propagation  Laboratory 
Boulder,  Colorado  80302 


1.  Introduction 

The  most  hazardous  aspects  of  aviation  are  the  landing  and  take- 
off when  both  airspeed  and  altitude  are  necessarily  small.  At  these 
times  aircraft  are  the  most  vulnerable  to  small  scale  meteorological 
phenomena.  In  recent  years,  the  hazards  of  low  level  wind  shear  have 
become  apparent  and  some  crashes  that  previously  would  have  been  attri- 
buted to  "pilot  error,"  are  now  recognized  as  due  to  wind  shear.  How 
many  crashes  are  wind  shear  related  is  not  known.  However,  Fujita  and 
Caracena  (1977)  analyzed  three  commercial  airline  crashes  that  occurred 
in  1975  and  1976  and  determined  that  they  were  due  to  wind  shear. 

There  are  several  remote  sensing  techniques  that  could  be  used  to 
measure  wind  shear  at  an  airport.  Which  one  is  "best"  depends  on  the 
requirements  for  an  operational  wind  shear  detection  system.  These 
requirements  depend  in  turn  on  the  structure  of  the  meteorological 
phenomenon  that  causes  the  hazardous  wind  shear.  Fujita  and  Caracena 
(1977)  found  the  hazardous  wind  shear  was  caused  by  an  intense,  localized 
downdraft,  which  they  term  a "downburst."  They  found  that  the  hazardous 
area  may  be  only  a few  km  across  and  can  move  horizontally  very  fast.  On 
the  basis  of  this,  a reasonable  requirement  for  a wind  shear  detection 
system  is  that  it  be  capable  of  measuring  shear  in  the  air  volume  in 
which  the  aircraft  will  be  vulnerable  to  wind  shear.  This  capability 
would  require  measurement  along  the  glide  path  from  the  mid-point  of  the 
runway  to  a range  where  the  aircraft  would  have  sufficient  altitude  to 
recover  from  a "downburst."  From  the  data  on  actual  crashes  and  near 


crashes  given  by  Fujita  and  Caracena  (1977),  a range  of  10,000  feet 
would  be  required.  An  adequate  safety  margin  might  increase  this  to 
15,000  ft  or  about  5 km. 

One  remote  sensor  that  has  potential  for  meeting  these  require- 
ments is  the  FM-CW  Doppler  radar.  (An  explanation  of  the  principle  behind 
this  type  of  radar  is  given  in  Appendix  A.)  The  Wave  Propagation  Laboratory 
FM-CW  radar  is  contained  on  two  trailers  which  can  be  relocated  in  one 
day.  The  system  uses  two  8-fooL  antennas  which  can  be  scanned  in  a 
hemisphere.  The  wavelength  is  10  cm  and  the  transmitted  power  is  200  W. 

When  hydrometeors  are  present,  they  act  as  tracers  of  the  wind;  how- 
ever, the  radar  has  sufficient  sensitivity  to  detect  backscatter  from 
the  refractive- index  fluctuations  in  the  clear  air.  These  turbulence- 
induced  fluctuations  then  act  as  tracers  of  the  wind.  Since  attenuation 
is  small  at  these  frequencies,  the  system  has  all-weather  capability. 

This  report  describes  an  experiment  designed  to  test  the  feasibility 
of  using  an  FM-CW  Doppler  radar  to  detect  hazardous  wind  shear.  The 
main  output  of  the  experiment  is  the  collection  of  strength-of-return 
data  as  a function  of  time  of  day  and  month,  and  these  data  are  presented  in 
such  a way  as  to  be  applicable  for  any  electromagnetic  backscatter  radar. 

An  explanation  of  how  the  radar  was  calibrated  to  obtain  these  data  is 
given  in  Section  2.  The  strength-of-return  data  are  presented  in  histogram 
form  in  Section  3 and  Appendix  D.  A second  goal  of  the  experiment  was 
to  attempt  radar  operation  at  low  elevation  angles  and  this  is  discussed 
in  Section  4.  The  results  of  three  other  related  activities  are  discussed 
in  Section  5 and  conclusions  are  presented  in  Section  6. 

2.  Radar  Calibration  and  Calculations 

It  is  generally  not  necessary  to  calibrate  an  FM-CW  Doppler  radar 
to  make  velocity  or  wind  shear  measurements.  However  calibration  is 
required  to  make  the  strength-of-return  measurements  reported  here. 

For  a dual  antenna  radar,  two  types  of  calibration  are  required.  First, 
the  electronic  portion  of  the  radar  must  be  calibrated  so  that  the  area 


2 


under  a Doppler  spectrum  can  be  related  to  the  received  power  associated 

with  that  spectrum.  Second,  since  the  effective  antenna  pattern  Is 

actually  the  overlap  between  the  transmitting  and  receiving  antennas, 

this  overlap  must  be  quantified.  The  next  two  sections  discuss  these 

calibrations  and  the  last  section  discusses  how  the  calibrated  return 

2 

power  is  related  to  C^,  the  atmospheric  parameter  of  interest  in  clear 
air  wind  measurement. 

2.1  Electronic  Calibration 

The  basic  idea  in  calibration  of  the  electronic  portion  of  our 
radar  is  to  inject  a known  test  signal  into  the  receiver.  Then  the 
proportionality  constant  between  received  power  and  radar  output  can 
be  determined.  Two  different  configurations  were  used  to  obtain  the  test 
signal.  The  first  was  to  take  the  test  signal  from  the  output  of  the 
high  power  transmitter  and  return  it  to  the  receiver  trailer  by  coaxial 
cable.  The  major  difficulty  with  this  arrangement  was  that  the  attenuation 
of  this  coaxial  cable  changed  with  temperature.  The  temperature  changes 
were  caused  mainly  by  the  sun  heating  the  cable.  To  avoid  this  problem 
we  took  the  test  signal  from  inside  the  trailer  so  that  any  coaxial 
cable  was  not  in  the  sun.  This  necessitated  using  a power  monitor  to 
measure  and  record  the  transmitted  power.  A block  diagram  of  the  cali- 
bration arrangement  we  used  most  of  the  time  is  shown  in  Fig.  2.1. 

The  calibration  procedure  (which  was  performed  usually  once  per  day) 
is  to  turn  off  the  TWT,  connect  the  test  signal,  and  then  record  one  record 
out  of  the  signal  processor.  Then  the  test  signal  is  disconnected  and 
a second  record  is  recorded  to  give  a noise  baseline  with  no  input.  One 
record  is  then  subtracted  from  the  other  and  the  difference  (in  the  10  us 
range  cell)  is  the  spectrum  of  the  test  signal.  Since  the  strength  of 
the  test  signal  is  known  (from  the  power  meter  reading)  and  the  area 
under  the  test  signal  spectrum  is  known,  the  proportionality  constant 
relating  spectral  area  to  received  power  can  be  found. 
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Figure  2.1.  FM-CW  radar  calibration  subsystem. 

The  only  difficulty  with  this  approach  Is  that  the  RF  low  noise 
amplifier  is  not  in  the  calibration  loop.  Gain  changes  in  this  amplifier 
were  kept  to  a minimum  by  a constant  temperature  heater.  The  gain  was 
checked  periodically  to  monitor  any  ’ong  term  changes.  The  main  cali- 
bration error  is  probably  due  to  unknown  gain  changes  in  this  RF  amplifier 
and  we  believe  that  the  unknown  changes  are  2 dB. 

2.2  Antenna  Calibration 

Since  the  FM-CW  radar  uses  seperate  antennas  for  transmitting  and 
receiving,  an  unconventional  calibration  and  measurement  procedure  must 
be  used.  First,  we  measured  antenna  patterns  using  conventional  antenna 
range  techniques.  The  patterns  for  the  two  antennas  were  very  nearly 
identical  and  one  of  the  measured  patterns  is  shown  in  Fig.  2.2.  Then 
we  used  an  analytical  argument  to  show  how  thi?  single  pattern  is  changed 
into  an  effective  antenna  pattern  that  is  a function  of  range  as  well 
as  angle. 
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Figure  2. 2.  Antenna  pattern  for  receiving  antenna. 
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The  antennas  beams  point  in  parallel  directions,  so  there  is  only 
a small  overlap  at  close  ranges  and  complete  overlap  at  long  ranges. 

■ 

This  overlap  is  characterized  by  multiplying  one  antenna  pattern  by  the 
other  and  integrating  the  result.  Assume  the  two  antennas  are  separated 
by  a distance  D and  are  pointed  in  parallel  directions  as  shown  in  Fig.  2.3. 

Define  F(R),  the  overlap  at  range  R,  as  the  integral  of  the  product  of 
the  two  antenna  patterns.  Using  small-angle  approximations  in  Fig.  2.3, 
we  see  that 


F(R)  = § / G 


■mm 


dW 


where  the  term  outside  the  integral  is  necessary  to  insure  the  boundary 
conditions  of  the  overlap.  We  require  that  F(R-*0)  = 0 and  FCR^00)  =1. 
Evaluating  the  constant  and  changing  the  variable  of  integration  we  obtain 


F(R) 


/Cl (9  + |g)G2(e  - S^)d8 
fG^d)  G2( 8)  d0 


But  this  is  just  the  normalized  cross-correlation  function  between  the 
two  antenna  patterns,  and  in  our  case,  since  the  two  antenna  patterns  are 
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Figure  2.  3.  Beam  overlap  for  tuo  antennas. 

equal,  this  becomes  the  autocorrelation  function,  R (•)•  So,  the  overlap 

G 

is  given  by 

F(R)  = Rg(|)  • 

The  overlap  can  be  calculated  with  lagged  products  and  is  shown  in  Fig.  2.4. 

In  any  return  strength  of  signal  calculation,  the  effects  of  the 
antenna  overlap  can  be  taken  into  account  by  simply  multiplying  either  the 
antenna  pattern  or  the  effective  area  by  F(R).  This  was  done  in  all  of 
our  calculations. 

2 

2.3  C Calculations 
n 

Once  the  radar  is  calibrated,  the  return  power  for  each  range  cell 
can  be  determined,  but  these  results  would  apply  only  to  our  radar.  A 
better  approach  is  to  calculate  a pertinent  atmospheric  parameter  that 
would  be  applicable  to  other  radars  as  well.  For  clear  air  radar  scattering 


6 


Figure  2.  4.  Range  antenna  function. 

the  best  atmospheric  parameter  is  the  structure  constant  of  radio  re- 

2 -2/3 

fractive- index,  which  is  generally  denoted  as  and  has  units  of  m .If 

the  atmospheric  turbulence  is  homogeneous  and  isotropic  and  n(rQ)  is 

2 

the  refractive- index  at  position  r , then  C is  defined  by 

o n 

[n(rQ)  - n(rQ  + r)]2  * C2  r2^3  , 

where  r is  the  separation  distance  between  two  points.  Whenever  the 
scale  size  of  interest,  X/2,  lies  within  the  inertial  subrange  of 
turbulence,  the  radar  cross-section  per  unit  volume,  n,  is  given  by, 
(Ottersten,  1969) 


n 


0.38  C2  X"1/3 
n 


The  inertial  subrange  lies  between  the  inner  and  outer  scale  of  turbulence, 

and  in  the  boundary  layer  the  inner  scale  is  on  the  order  of  a mm.  The 

outer  scale  is  determined  by  atmospheric  stability;  for  high  stability 

this  outer  scale  is  about  5 m and  under  convective  conditions,  it  increases 

2 

to  perhaps  2 km.  So  the  C values  given  here  should  be  useful  for  any 

n 

conceivable  radar  application  in  the  boundary  layer. 
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To  derive  the  relation  between  C and  radar  output,  we  start  with 

n 

the  conventional  radar  equation  for  distributed  targets.  Let  P^  be 
the  power  at  the  output  of  the  receiving  antenna;  PT  be  the  transmitted 
power  at  the  input  to  the  transmitting  antenna;  be  the  effective  area 
of  one  of  the  identical  antennas;  A be  the  range  cell  size;  r|  be  the 
radar  cross  section  per  unit  volume;  and  R be  the  range  of  interest. 

Then  we  have 


P 

r 


p_  a An 

T e 


4tt  R 


However,  the  overlap  of  the  antenna  beams  discussed  in  the  previous  section 
must  be  included  to  adjust  for  lack  of  intersection  at  lower  ranges.  This 
gives 


P 

r 


PT  RG  (f>  Ae  a 1 
4tt  R2 


Solving  for  after  substituting  for  n we  obtain 


c2.!i 

" PT 


33 


R2  A1/3 


Vr> 


Now,  since  the  target  fills  the  beam,  we  must  make  a correction  since 
the  antenna  gain  at  the  edges  of  the  beam  is  different  than  that  at  the 
center  (Probert- Jones,  1962).  Also  we  use  X = .1  m and  A^  = 2.52  m to 
obtain 


8.43  Pr 


?T 


R_ 

A 


The  only  thing  left  is  to  insert  the  calibration  of  the  radar  receiver 
and  signal  processor  so  that  P can  be  replaced  by  a quantity  from  the 
output  of  the  radar.  Let  A^  be  that  quantity  which  is  equal  to  the  area 
under  the  velocity  spectrum  for  the  range  of  interest.  Also  let  Ag  be 
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the  area  under  the  velocity  spectrum  of  the  test  signal  which  has  power 

P . It  is  then  clear  that 
s 


and  the  equation  for  C becomes 

n 


c2  = 

n 


8.43 

Vf> 


P ? 

s R , _ .-2/3 

p-  (meters) 

T 


Note  that  the  units  on  A and  A are  unimportant  since  they  cancel.  This 

r s 2 

is  the  equation  that  was  used  to  calculate  the  measured  values  given 
in  the  next  section. 


3.  C Measurements 
n 

2 

The  radar  made  measurements  from  March,  1977  through  February 
1978.  During  the  spring  and  summer,  we  attempted  to  take  data  continuously, 
and  we  failed  only  during  periods  of  transporting  the  radar  or  during 
periods  of  equipment  failure.  During  the  late  fall  and  winter,  we  took 
data  on  an  intermittent  basis  because  the  radar  computer  was  being  used 
for  processing  data. 


3.1  Mode  of  Operation 

During  the  entire  experiment,  we  deviated  only  slightly  from  a fairly 
standard  mode  of  operation.  The  radar  antennas  were  pointed  at  an  elevation 
angle  of  60o-75°  and  the  azimuth  was  any  convenient  value,  since  the  height 
and  height  resolution  is  unaffected  by  the  azimuth  angle.  The  frequency 
was  generally  just  a few  MHz  below  3 GHz,  and  the  transmitted  power  was 
approximately  200  W.  The  sweep  length  was  generally  3 ms,  and  50  sweeps 
were  coherently  processed.  From  Appendix  A we  see  that  this  means  the 
output  has  a 50  point  velocity  spectrum  for  each  of  10  range  cells. 
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Because  of  power  line  harmonics,  we  did  not  attempt  to  process  the  data 

in  the  lowest  range  bin  and  this  left  9 range  cells  of  length  161  m each. 

The  spectra  were  recorded  on  magnetic  tape,  and  the  calculations  to  extract 
2 

values  were  performed  off-line. 

A record  of  data  consists  of  500  data  points  (ten  50-point  velocity 

spectra)  and  one  record  was  recorded  every  60  seconds.  Each  record  was 

a running  average  of  the  radar  output  over  approximately  60  seconds.  In 

2 

the  off-line  processing  programs  the  value  for  each  range  cell  was 

calculated  for  each  of  the  1-minute  records  as  described  in  the  previous 

2 

section.  This  resulted  in  60  C measurements  at  that  range  for  each 

n 

hour  of  operation,  and  these  60  measurements  were  used  to  form  a histogram 
2 

of  values.  The  histograms  were  written  onto  a compressed  data  tape 
which  was  used  to  study  histograms  as  functions  of  height  or  time  or 
hour  of  the  day. 

3.2  Time  Ordered  Histograms 

In  the  Boulder,  Colorado,  area  over  the  past  two  years  we  have 
2 

observed  that  C in  the  boundary  layer  can  vary  over  a 63  dB  range. 

n 2 

Because  of  this  large  range,  the  histogram  variable  is  -log  0 C rather 

2 2 iU  n 
than  C . For  each  decade  of  C values,  three  histogram  bins  were 

n 2 n -15  -15  -2/3 

formed.  For  example,  a C bin  extends  from  1 x 10  to  2.15  x 10  m 

n -15  -15  -2/3 

The  next  higher  bin  extends  from  2.15  x 10  to  4.64  x 10  m and  the 

-15  -15  -2/3 

next  higher  from  4.64  x 10  to  10  x 10  ra  . This  results  in  three 

equal  sized  bins  on  a logarithmic  scale  and  is  repeated  for  each  decade 
2 

of  C values, 
n 

Figure  3.1  is  a sequence  of  hourly  histograms  that  extend  over  r 

two-week  period.  The  five  columns  of  numbers  on  the  left  are  month. 

day,  year,  hour,  and  minute.  The  row  of  numbers  extending  from  11  to 
2 

19  are  -log,-  C where  each  decade  has  been  broken  into  three  bins  as 
°10  n 

described  above.  The  1-hour  histograms  are  shown  as  rows  across  the 

page  where  each  number  represents  how  many  1-minute  measurements  of 
2 

C fell  into  that  bin  of  values  during  the  hour.  This  figure  shows  the 

n 2 
hour-to-hour  and  day-to-day  variation  of  a typical  set  of  C data. 


10 


RANCE  FOR  THIS  CELL  IS  = 805  METERS  ***«« 

DATE  MONTH/DAY/YEAR  5/14/77 
TIME  HOURS  MINUTES  SECONDS  11  41  55 


10  0 

4 64  2 

1 15 

FILE 

9 

REC 

38 

11  11 

11  12  12  12  13  13 

13 

14 

I 4 

14 

15 

15 

15 

16  16 

16 

17 

5 

14 

77 

11 

41 

7 

19 

15 

16 

3 

5 

14 

77 

12 

42 

10 

13 

23 

12 

1 

1 

5 

14 

77 

13 

42 

44 

16 

5 

14 

77 

14 

43 

2 

1 1 

8 

6 

30 

3 

5 

14 

77 

15 

43 

10 

13 

37 

5 

14 

77 

16 

43 

15 

9 

22 

14 

5 

14 

77 

17 

44 

1 

9 

7 

8 

20 

11 

4 

5 

14 

77 

18 

44 

2 

37 

21 

5 

14 

77 

19 

45 

40 

19 

1 

5 

14 

77 

20 

47 

19 

29 

9 

3 

5 

14 

77 

21 

49 

1 

5 

7 

23 

20 

3 

1 

5 

14 

77 

22 

51 

4 

26 

14 

13 

2 

1 

5 

14 

77 

23 

53 

7 

9 

18 

8 

12 

6 

5 

15 

77 

0 

54 

12 

37 

10 

1 

5 

15 

77 

1 

55 

20 

37 

3 

5 

15 

77 

2 

57 

3 

15 

6 

16 

9 

6 

4 

5 

15 

77 

4 

0 

19 

17 

5 

4 

6 

6 

2 

1 

5 

15 

77 

5 

2 

5 

41 

10 

3 

1 

5 

15 

77 

6 

4 

1 

33 

20 

6 

5 

15 

77 

7 

6 

16 

34 

9 

1 

5 

15 

77 

8 

9 

1 1 

45 

4 

5 

15 

77 

9 

10 

23 

26 

2 

5 

4 

5 

15 

77 

10 

10 

4 

37 

7 

4 

2 

2 

1 

1 

1 

5 

15 

77 

11 

10 

3 

38 

7 

2 

7 

2 

1 

5 

15 

77 

12 

10 

41 

15 

1 

1 

2 

5 

15 

77 

13 

1 1 

1 

50 

9 

5 

15 

77 

14 

11 

14 

35 

10 

1 

5 

15 

77 

15 

11 

13 

31 

15 

1 

5 

15 

77 

16 

12 

17 

23 

17 

3 

5 

15 

77 

17 

12 

4 

15 

40 

1 

5 

15 

77 

18 

12 

1 1 

46 

3 

5 

15 

77 

19 

13 

3 

47 

10 

5 

15 

77 

20 

13 

1 

10 

4 

5 

33 

7 

5 

15 

77 

21 

13 

23 

12 

24 

1 

5 

15 

77 

22 

14 

4 

48 

8 

5 

15 

77 

23 

14 

8 

39 

13 

5 

16 

77 

0 

14 

1 

19 

29 

1 1 

5 

16 

77 

1 

15 

28 

32 

5 

16 

77 

2 

15 

6 

40 

12 

1 

5 

16 

77 

3 

15 

58 

2 

5 

16 

77 

4 

16 

24 

34 

1 

1 

5 

16 

77 

5 

16 

1 

26 

21 

9 

3 

5 

16 

77 

6 

18 

3 

28 

27 

2 

5 

16 

77 

7 

19 

3 

57 

5 

16 

77 

8 

20 

35 

25 

5 

16 

77 

9 

20 

2 

51 

6 

1 

5 

16 

77 

10 

20 

1 

43 

16 

5 

16 

77 

11 

21 

23 

36 

1 

5 

16 

77 

12 

21 

19 

38 

2 

1 

5 

16 

77 

13 

21 

36 

20 

3 

5 

16 

77 

14 

22 

32 

28 

5 

16 

77 

15 

22 

34 

24 

2 

5 

16 

77 

16 

22 

4 

29 

27 

5 

16 

77 

17 

23 

27 

33 

5 

16 

77 

18 

23 

1 

48 

11 

. 2 
Figure  3.1(a).  Hourly  histograms  of  C at  805  m. 


11 


T 


DATE 

TIME 


» RANGE  FOR  THIS  CELL  IS  » 
MONTH/DAY/YEAR  5/16/77 
HOURS  MINUTES  SECONDS  19  23 
10  0 4 64 

11  11  11  121212 


805  METERS 


53 


2 15 


5 

16 

77 

19 

23 

5 

16 

77 

20 

24 

5 

16 

77 

21 

24 

5 

16 

77 

22 

24 

5 

16 

77 

23 

25 

5 

17 

77 

0 

25 

5 

17 

77 

1 

26 

5 

17 

77 

2 

26 

5 

17 

77 

3 

27 

5 

17 

77 

4 

29 

5 

17 

77 

5 

31 

5 

17 

77 

6 

33 

5 

17 

77 

7 

34 

5 

17 

77 

8 

35 

5 

17 

77 

9 

35 

5 

17 

77 

10 

35 

5 

17 

77 

1 1 

36 

5 

17 

77 

12 

36 

5 

17 

77 

13 

16 

5 

17 

77 

14 

37 

5 

17 

77 

15 

37 

5 

17 

77 

16 

37 

5 

17 

77 

17 

38 

5 

17 

77 

18 

38 

5 

17 

77 

19 

38 

5 

17 

77 

20 

39 

5 

17 

77 

21 

39 

5 

17 

77 

22 

39 

5 

17 

77 

23 

40 

5 

18 

77 

0 

40 

5 

18 

77 

1 

40 

5 

18 

77 

2 

41 

5 

18 

77 

3 

41 

5 

18 

77 

4 

41 

5 

18 

77 

5 

43 

5 

18 

77 

6 

45 

5 

18 

77 

7 

46 

5 

18 

77 

8 

46 

5 

18 

77 

9 

46 

5 

18 

77 

10 

49 

5 

18 

77 

1 1 

49 

5 

18 

77 

12 

49 

5 

18 

77 

13 

50 

5 

18 

77 

14 

50 

5 

18 

77 

15 

50 

5 

18 

77 

16 

50 

5 

18 

77 

17 

51 

5 

18 

77 

18 

51 

5 

18 

77 

19 

51 

5 

18 

77 

20 

52 

5 

18 

77 

21 

52 

5 

18 

77 

22 

52 

5 

18 

77 

23 

53 

5 

19 

77 

0 

53 

5 

19 

77 

1 

53 

52 

47 


FILE 

9 

14 

14 

14 

15 

5 

53 

2 

42 

7 

35 

17 

3 

3 

7 

42 

6 

15 

40 

5 

3 

29 

27 

1 

4 

22 

28 

2 

23 

34 

2 

20 

36 

4 

1 1 

32 

15 

3 

9 

29 

18 

2 

58 

30 

30 

1 

9 

46 

4 

14 

45 

1 

2 

56 

2 

4 

56 

25 

34 

1 

32 

28 

54 

5 

34 

2 

1 

9 

11 

29 

5 

34 

26 

43 

17 

12 

37 

6 

2 

4 

14 

24 

12 

4 

14 

33 

7 

11 

26 

9 

8 

34 

7 

1 

25 

4 

6 

3 

21 

2 

42 

3 

5 

30 

24 

1 

17 

22 

21 

29 

5 

15 

5 

39 

16 

26 

24 

8 

21 

3 

24 

19 

36 

7 

27 

9 

29 

15 

34 

6 

15 

4 

3 

6 

10 

3 

21 

5 

37 

11 

14 

41 

2 

30 

1 

9 

3 

16 

24 

REC 


1 1 


94 


17 


2 22  25  10 
12  39  6 

18  15  8 1 


35 

4 

4 

26 

39 


Figure  3.1(b).  Hourly  histograms  of  at  80S  m. 
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Figure  3.1(c).  Hourly  histograms  of  C at  805  m. 
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Figure  3.1(e).  Hourly  histograms  of  C at  805  m. 
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Figure  3.1(f).  Hourly  histograms  of  C at  805  m. 
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3. 3 Combined  Monthly  Histograms 


The  monthly  variation  of  can  be  found  by  combining  hourly 

histograms  for  each  month  and  these  are  shown  in  Appendix  D.  These 

histograms  contain  the  effects  of  birds,  insects,  airplanes,  rain,  and 

snow  storms,  etc.  Thus,  the  upper  (leftmost)  part  of  the  histogram 

can  not  be  attributed  completely  to  clear  air  return.  However,  it  is 

likely  that  the  lower  part  of  the  histogram  is  relatively  free  of  these 

2 

other  targets.  Figure  3.2  shows  C^  data  obtained  by  integrating  the 

lower  part  of  the  histograms  to  obtain  three  values,  one  that  will  be 
2 

greater  than  C 5%  of  the  time,  a second  that  will  be  greater  than 

2 n 2 
C 30%  of  the  time,  and  a third  that  will  be  greater  than  C 50%  of 
n n 

the  time.  In  Figure  3.2  these  levels  are  plotted,  for  each  of  the  nine 
range  cells,  as  a function  of  month. 

One  surprising  feature  of  Figure  3.2  is  the  drop  from  June  to  July, 
since  one  would  expect  similar  returns  in  June  and  July.  One  possible 
explanation  for  this  drop  is  that  the  radar  was  moved  to  a different 
location.  The  March  to  June  data  were  taken  near  Boulder,  Colorado, 
near  large  areas  of  irrigated  farmland.  The  July  to  October  data  were 
taken  at  Buckley  Field  east  of  Denver.  There  is  little  or  no  irrigation 
nearby  and  the  ground  cover  is  mostly  high  plains  grass  and  sagebrush.  This 

lack  of  ground  level  moisture  may  account  for  the  June  to  July  drop  in 

2 2 
measured  C . If  this  is  true,  it  is  likely  that  values  of  C in  more 
n J n 

humid  regions,  such  as  the  central  and  eastern  United  States,  are  consider- 
ably greater  than  the  measurements  presented  here. 

A second  surprising  feature  of  Figure  3.2  is  the  difference  between 
the  measurements  of  March  1977  and  February  1978.  Normally,  one  would 
not  expect  them  to  be  greatly  different.  However  the  winters  of  76-77 
and  77-78  were  drastically  different.  In  Colorado,  the  winter  of  76-77 

was  unusually  mild  whereas  77-78  was  more  normal  and  had  much  colder 

2 

temperatures.  The  warmer  winter  resulted  in  higher  C^  values  than 
normal  and  this  is  shown  in  our  measurements. 
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Figure  2.2(a).  Monthly  C data  at  161  m.  Graphs  show  C value  that 

^2  ** 
was  greater  than  observed  C ‘ for  that  percentage  of  time. 
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Figure  3.2(f).  Monthly  C ^ data  at  966  m. 


TABLE  3.1.  Time  Operated  per  Month 


t 


l 


i 


[' 


Month 

Minutes 

Days 

March 

19260 

20 

April 

26640 

23 

May 

28140 

22 

June 

17640 

21 

July 

12060 

11 

August 

14220 

17 

September 

26160 

22 

October 

3660 

5 

November 

6540 

5 

December 

7740 

7 

January 

3960 

4 

February 

3600 

4 

As  one  would  expect,  is  significantly  lower  in  the  winter 

months.  However,  this  should  have  little  effect  on  an  operational  wind 
shear  system,  because  hazardous  wind  shear  generally  does  not  occur  in 
the  winter. 

The  number  of  minutes  operated  each  month  is  shown  in  Table  3.1. 
The  total  operating  time  was  2827  hr  and  covered  160  days. 

3.4  Combined  Hour-of-the-Day  Histograms 

The  data  can  also  be  viewed  on  an  hour-of-the-day  basis  for  each 
range  cell,  but  the  results  for  just  the  805  m range  cell  will  be  pre- 
sented here.  The  histograms  for  each  of  the  24  hours  in  the  day  are 

shown  in  Appendix  D.  The  5%,  30%  and  50%  levels  of  these  histograms 

2 

are  shown  in  Figure  3.3  which  has  two  peaks  in  the  values.  The 
first  peak  occurs  about  midnight  and  the  second  at  midday.  These  peaks 
are  probably  due  to  a combined  effect  of  the  temperature  gradient  and 
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Figure  3.3.  Hourly  data  at  805  m.  Graphs  show  value  that 
was  greater  than  observed  cj1  for  that  percentage  of  time. 


mechanical  turbulence.  At  sunrise  there  are  large  gradients,  but  very 

little  turbulence.  As  the  sun  rises  higher,  solar  heating  and  the 

2 

resulting  convection  cause  increased  which  reaches  a peak  shortly 

after  noon.  There  is  considerable  convection,  but  very  little  departure 

from  an  adiabatic  temperature  gradient.  Then  the  heating  is  decreased 

2 

in  the  afternoon  and  the  values  decrease.  As  the  ground  cools  in 

late  afternoon  and  evening,  the  atmospheric  stability  increases  and  this 

causes  large  gradients.  These  large  gradients  require  only  small  amounts 

2 2 

of  turbulence  to  cause  large  C values,  and  C reaches  a peak  about 

n n 

midnight.  After  midnight,  the  stability  increases,  but  the  turbulence 
2 

dies  out  so  that  C decreases  to  a minimum  just  before  sunrise  when  the 
n J 

cycle  starts  over  again. 

So,  the  two  peaks  in  Figure  3.3  are  due  to  two  different  effects. 

The  midnight  peak  is  due  to  high  stability,  the  resulting  large  gradients, 

and  a small  amount  of  turbulence.  The  noon  peak  is  due  to  strong  con- 

2 

vection  and  small  gradients.  The  minimum  in  C at  sunrise  is  due  to  a 

n 

stable  atmosphere  (with  large  gradients)  but  very  little  turbulence. 


3. 5 Probability  Law  for  C 


An  item  of  importance  for  system  design  is  a general  probability 

2 2 
law  for  C . Because  large  and  small  values  of  C occur  with  some 
n 2 n 

regularity,  it  is  likely  that  does  not  follow  a Gaussian  distribution. 

A more  suitable  model  for  a parameter  with  a large  variation  is  the 

2 

lognormal  distribution,  i.e.,  log  C being  Gaussian.  A straightforward 

n j 

way  to  check  this  is  to  plot  log  C on  probability  paper.  If  the  plot 

2 n 

is  a straight  line,  then  C follows  a lognormal  distribution. 

n 2 

Figure  3.4  shows  plots  of  log  C at  805  m for  each  month  data  were 

taken.  Figure  3.5  is  a similar  plot  except  the  variable  plotted  is 
2 

log  C for  every  third  hour  of  the  diurnal  cycle.  Only  the  March  to 
n 

October  data  were  used  for  Fig.  3.5.  The  conclusion  from  these  two 
2 

figures  is  that  can  be  assumed  to  be  lognormal.  Further,  since  the 
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Figure  3.4.  Probability  distribution  of  monthly  log  Cn  values. 
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Percentage  of  Time  That  Cp  Is  Less  Than  Abcissa 
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Percentage  of  Time  That  Cp  Is  Less  Than  Abcissa 


slopes  of  the  lines  in  these  two  figures  are  fairly  constant,  the  standard 
2 

deviation  of  log  C can  be  taken  as  a constant  value  of  about  6 dB. 

n 2 
This  is  important  because  it  simplifies  the  lognormal  model  of  from 

a two  parameter  model  to  a one  parameter  model,  i.e.,  only  the  mean  of 
2 

log  needs  to  be  specified. 

The  specific  value  one  should  use  for  the  mean  of  log  C ^ depends 

n 

on  the  month,  hour  of  day,  height,  as  well  as  location.  The  data  here, 

i.e.,  the  50%  lines,  could  be  used  for  a high  plains  meteorological  regime. 

2 

Oossard  (1976)  presents  a method  for  transferring  profiles  of  data  from 

one  air  mass  to  another.  Using  this  technique,  one  could  get  a rough 

2 

idea  of  the  mean  of  log  at  other  locations. 

3.6  Comparison  with  in-situ  Measurements 

We  have  compared  our  measurements  with  three  much  smaller  data  sets 

taken  under  similar  conditions,  and  found  reasonable  agreement.  Our 

data  set  is  from  the  boundary  layer,  but  above  the  surface  layer.  There 

are  numerous  data  sets  in  the  surface  where  models  of  height  dependence 

2 

exist.  However,  we  encountered  difficulty  taking  calibrated  mea- 
surements in  the  surface  layer  so  we  can  not  compare  with  surface  layer 
measurements.  (There  is  no  difficulty  making  wind  measurements  in  the 
surface  layer  with  our  radar  since  a power  calibration  is  not  required.) 

Of  the  three  data  sets,  we  have  compared  with,  two  use  radio  refrac- 
tive index  as  we  do  and  the  third  uses  only  the  temperature  fluctuations. 
However,  this  third  one  was  taken  during  a relatively  dry  month  so  that 

hopefully  temperature  is  the  contributing  factor  to  the  radio  refractive- 

2 

index.  Gossard  (1976)  shows  the  relative  contributions  to  C for  tempera- 

n 

ture  and  humidity. 

The  first  data  set  (Bean  et  al.,  1970)  is  four  tower  measurements 
of  refractivity  fluctuations  using  two  microwave  refractometers  spaced 
20  cm  apart.  The  measurements  were  made  in  southeastern  Colorado  in 
October  1970.  The  height  range  is  75-130  m and  the  four  measurements 
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range  from  1.2  x 10  ^ to  8.5  x 10  ^ m The  measurements  were  taken 

when  the  ref ractometers  were  In  regions  with  strong  temperature  gradients 

and  so  it  is  not  surprising  that  they  are  5 to  23  dB  higher  than  our 

mean  value  for  October  at  161  m.  Figure  2 of  Bean  et  al.  (1970)  shows 
2 

that  the  peak  value  in  the  region  of  strong  temperature  gradients 
is  about  15  dB  above  a baseline.  The  conclusion  here  is  that  the 
agreement  is  reasonable. 

2 

The  second  data  set  (Gossard,  1976)  is  estimates  of  C based  on 

n 

rawinsonde  measurements  and  temperature  fluctuation  models.  The  relevent 

figure  is  17  where  the  summer  values  are  seen  to  be  somewhat  greater  than 

-15  -2/3  -15  -2/3 

10  ra  , while  the  winter  values  are  somewhat  less  than  10  m 

Our  values  are  less  than  this  by  3 to  7 dB  and  since  Gossard 's  data  set 

is  for  Omaha,  Oklahoma  City,  and  Boston,  the  difference  is  not  surprising. 

The  third  data  set  (Ochs  and  Lawrence,  1972)  is  from  temperature 

measurements  made  by  aircraft  in  southeastern  Colorado  during  November 

1971.  If  the  humidity  was  quite  low,  these  measurements  should  not  be 

greatly  different  than  ours.  The  relevent  figure  is  the  average  profile 

2 

over  land  on  page  36,  where  the  C range  for  heights  above  160  m,  is  from 

_17  —16  ^ 

5 x 10  to  5 x 10  . This  is  in  good  agreement  with  our  values  for 

the  month  of  November. 

The  agreement  between  our  data  set  and  these  three  others  gives  some 
confidence  that  our  data  is  not  drastically  in  error. 

4.  Operation  at  Low  Elevation  Angles 

If  any  remote  sensing  shear  detection  system  is  to  be  effective, 
it  must  detect  hazardous  shear  near  the  region  where  the  aircraft  is 
the  most  vulnerable,  i.e.,  the  glidepath.  The  spatial  structure  of 
meteorological  phenomena  that  resulted  in  hazardous  wind  shear  on  two 
occasions  has  been  studied  by  Fuj ita  and  Caracena  (1977).  They  found 
that  the  hazardous  region  can  be  as  small  as  a few  km  wide  and  that  this 
hazardous  region  moves  horizontally  very  fast.  Thus,  it  is  clear  that  a 
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shear  detection  radar  must  be  capable  of  continually  monitoring  the  wind 
at  a fairly  low  elevation  angle.  Further,  the  radar  must  be  capable  of 
measuring  winds  at  close  ranges  because  the  hazardous  shear  region  may 
be  positioned  over  the  runway.  This  combination  of  operation  at  low 
elevation  angles  and  close  ranges  is  difficult  for  most  radars.  A com- 
parison of  the  difficulty  for  a pulse  and  CW  radar  is  discussed  in 
Appendix  C.  This  section  will  discuss  these  difficulties  and  some 
possible  solutions  for  the  FM-CW  radar.  Also,  a way  of  displaying 
radial  shear  by  means  of  a Plan  Shear  Indicator  (PSI)  is  discussed  and 
examples  of  clear  air  PSI  displays  are  given.  It  is  clear  that  low 
elevation  angle  PSI  displays  from  the  center  of  the  runway  to  a range  of 
5 km  would  have  provided  adequate  information  to  avert  the  two  wind 
shear  caused  crashes  studied  by  Fuj ita  and  Caracena  (1977). 

4.1  Plan  Shear  Indicator  (PSI)  Display 

The  PSI  display  is  a qualitative  display  of  wind  shear  first  proposed 
by  Armstrong  and  Donaldson  (1969) . It  operates  like  a conventional  PPI 
in  that  one  can  readily  tell  the  azimuth  angle  and  range  to  an  area  of 
significant  shear.  The  velocity  spectra  for  increasing  ranges  are 
arranged  side-by-side  and  then  used  to  intensity  modulate  a PPI  oscillo- 
scope at  the  current  antenna  azimuth  angle.  Then  the  antenna  is  moved 
and  the  process  repeated.  If  the  only  return  is  from  stationary  targets, 
the  PSI  is  a set  of  arcs  of  concentric  circles  and  the  intensity  of  the 
arcs  is  proportional  to  the  received  power  from  that  range.  If  the 
targets  are  moving,  they  are  offset  from  the  arcs  and  the  offset  is 
proportional  to  the  velocity.  A change  in  velocity  from  one  range  cell 
to  the  next  (wind  shear  if  the  targets  are  clear  air)  is  then  quite 
apparent  because  the  spectra  do  not  have  the  same  spacing  as  the  ground 
clutter. 

A pulse-Doppler  radar  requires  special  signal  processing  to  produce 
a PSI;  however,  the  output  of  the  FM-CW  radar  is  already  in  the  required 
format  for  a PSI  without  additional  processing.  Appendix  A is  a simplified 
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explanation  of  FM-CW  Doppler  proc*-<  ing  an.’  It  hows  that  the  spectra 
occur  naturally  in  a range-ordered  i ' 1!.  f . lion.  So,  to  produce 

a PSI  display,  all  that  is  required  t . i z«  the  output  of  the 

spectrum  analyzer  to  a PPI  type  displa  ' <-  isplays  diown  here 

were  obtained  in  real  time  without  1 

Figure  A. 1(a)  shows  the  output  r . • trar.  malyzer  for  a 

summer  day  with  fairly  high  winds.  Th«  trr  w markers  show  tiie  location 
of  the  zero  velocity  points  (ZVP) , and  tin  wider  peaks  are  the  clear  air 
velocity  spectra.  The  range  cell  size  is  160  m,  the  maximum  unambiguous 
velocity  is  8.7  m/s,  and  the  elevation  angle  is  about  50°.  In  Fig.  4.1(a), 
the  closest  range  cell  is  to  the  left  and  the  farthest  is  to  the  right. 

An  incoming  velocity  displaces  to  the  left  of  the  ZVP  and  an  outgoing 
velocity  displaces  to  the  right.  Figure  4.1(b)  is  a PSI  display  taken 
under  the  same  conditions.  The  azimuth  scan  rate  is  about  l°/sec. 

It  is  immediately  obvious  that  the  wind  is  blowing  from  the  lower  left 
to  the  upper  right.  At  the  lower  part  of  the  PSI,  the  return  is  saturated, 
perhaps  because  of  an  aircraft.  Along  the  left  edge  of  the  PSI,  the  second 
and  third  range  cells  have  double  peaked  spectra.  The  double  peaked 
spectra  imply  two  regions  in  the  same  range  cell  with  different  velocities. 
This  would  imply  strong  wind  shear,  but  this  is  not  the  type  of  situation 
one  would  expect  to  be  hazardous.  However,  there  would  probably  be 
considerable  mechanical  turbulence  in  the  region  of  the  3rd  through  6th 
range  cells  along  the  left  side  of  the  region  covered  in  Fig.  4.1(b). 

The  maximum  azimuthal  scan  rate  for  our  radar  is  3°/sec,  and  when 
the  antenna  is  rotating  at  this  rate  the  averaging  time  must  be  reduced 
so  the  result  is  not  smeared.  Figure  4.2  shows  two  PSI  displays  taken  at 
a scan  rate  of  3° /sec  With  only  a 50  ms  averaging  time  as  compared  with 
about  3s  averaging  time  for  Fig.  4.1.  Again  the  elevation  is  about  50°. 
Figures  4.1  and  4.2  were  taken  on  different  days. 
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Figure  4.1(a).  A-seope  radar  output. 
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Figure  4.1(b).  Plan  Shear  Indicator  (PSI)  display. 
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Figure  4.3  is  a relatively  low  elevation  angle  PST  taken  alongside 
the  runway  at  Buckley  Field.  The  elevation  angle  of  29°  is  about  the 
minimum  we  could  achieve  without  removing  the  ground  clutter.  The  range 
cell  of  353  and  maximum  range  of  about  2800  m are  the  horizontal  distances, 
but  the  8 m/s  represents  a radial  velocity.  The  radial  lines  that  are 
white  all  the  way  out  are  due  to  signal  processor  saturation  by  air- 
craft or  ground  clutter.  At  a horizontal  range  of  2800  m and  an  elevation 

angle  of  29°,  the  height  being  probed  is  about  1550  m.  The  region  being 

probed  gets  higher  as  the  range  increases  and  eventually  is  above  the 
top  of  the  boundary  layer.  Because  of  the  drier  air,  the  return  is 
less  from  above  the  boundary  layer.  At  lower  elevation  angles  the  maximum 
range  should  increase.  So,  the  problem  with  achieving  ranges  of  5 km 
is  not  sensitivity,  but  rather  the  capability  to  operate  at  low  elevation 
angles. 

A PSI  type  of  display  can  give  a very  clear  picture  of  radial 

winds  along  the  glide  slope  if  the  problem  of  saturation  by  ground  clutter 

can  be  solved.  Tnis  is  covered  in  the  next  section. 

4.2  Ground  Clutter  Suppressors 

As  pointed  out  above,  operation  at  elevation  angles  near  the  glide 
slope  is  highly  desirable.  However,  the  ground  clutter  then  becomes 
large  and  saturates  the  signal  processor.  It  is  important  to  distinguish 
between  two  types  of  saturation,  receiver  saturation  and  signal  processor 
saturation.  Receiver  saturation  occurs  when  the  return  RF  signal  exceeds 
the  dynamic  range  of  either  the  RF  amplifier  or  the  mixer.  Whenever 
ground  clutter  causes  this  type  of  saturation,  wind  measurements  are 
not  possible.  The  approach  to  this  problem  would  be  to  use  antennas 
with  lower  sidelobes.  The  second  type  of  saturation  occurs  when  the 
ground  clutter  return  saturates  the  spectrum  analyzer.  Again,  wind 
measurements  are  not  possible  when  the  spectrum  analyzer  is  saturated. 
However,  this  large  return  is  due  to  ground  clutter  and  hence  repeats 
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Figure  4.3.  PSI  display  taken  at  Buckley  Field.  Elevation  angle  is  23° 
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from  one  sweep  to  the  next.  A properly  designed  ground  clutter  suppressor 
can  then  eliminate  the  periodic  component  and  greatly  reduce  the  dynamic 
range  of  the  signal  without  decreasing  the  wind  measuring  capability. 

We  have  designed,  built,  and  tested  a ground  clutter  suppressor  for 
the  FM-CW  radar.  On  the  basis  of  experiences  with  the  first  model,  we 
are  building  a second  unit  that  should  be  satisfactory.  The  first  model 
had  two  phases  of  operation.  During  the  acquisition  phase,  the  analog 
signal  was  digitized  and  these  samples  were  averaged  over  16  sweeps. 

This  sequence  of  stored  digital  samples  then  approximately  represented 
the  periodic  component.  During  the  cancellation  phase,  the  stored 
sequence  was  passed  through  a D-A  converter;  the  output  signal  was  smoothed 
and  then  subtracted  from  the  incoming  analog  signal.  The  result  of  the 
subtraction  is  free  of  a periodic  component.  The  advantage  of  this 
approach  is  that  the  analog  subtraction  introduces  no  quantizing  noise 
into  the  non-periodic  component.  However,  the  disadvantage  is  that  the 
unit  cannot  adapt  and  change  the  stored  signal  as  conditions  change. 

Also,  it  was  clear  that  more  flexibility  was  needed  in  the  number  of 
sweeps  averaged.  We  felt  a redesign  of  the  first  unit  was  needed,  rather 
than  modification. 

A block  diagram  of  the  redesigned  clutter  suppressor  is  shown  in 
Fig.  4.4.  The  basic  concept  is  similar  to  the  original  unit.  The  main 
difference  is  that  now  the  averaging  and  update  rate  is  controlled  by 
a micro-processor,  giving  the  added  flexibility  that  is  desirable  for  a 
research  device.  In  addition  to  block  averaging,  exponential  weighted 
averaging  can  be  used.  With  this  type  of  averaging,  the  stored  version 
of  the  periodic  component  is  continually  updated,  so  the  unit  actually 
adapts  to  changing  conditions. 

A general  criterion  for  clutter  suppressors  is  the  suppression 
range,  or  the  amount  the  clutter  will  be  reduced.  Since  the  output  D/A 
converter  has  8 bits,  the  maximum  that  clutter  can  be  suppressed  is  48 
dB  (6  dB/bit).  Whether  or  not  this  can  be  achieved  depends  on  the  periodicity 
of  the  clutter  signal.  The  experimental  work  with  the  first  clutter 
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Figure  4.4.  Ground  olutter  suppressor. 


suppressor  showed  clearly  that  the  clutter  signal  was  not  as  periodic  as 
had  been  assumed.  The  cause  of  the  non-periodicity  was  traced  to  an 
RF  instability,  the  elimination  of  which  is  discussed  in  the  next  section. 

4.3  Sweep  Frequency  Stabilizer 

An  instability  in  the  frequency  sweep  of  an  FM-CW  radar  lowers  the 
potential  clutter  suppression.  This  is  analyzed  in  Appendix  B,  and  from 
Fig.  B-2  the  frequency  stability  required  to  achieve  a given  clutter 
suppression  level  can  be  determined.  The  signal  source  used  in  the 
FM-CW  radar  is  an  yttrium- iron-garnet  (YIG)  transistor  oscillator  modulated 
by  a sweep  generator.  Noise  in  the  driving  stage  of  the  YIG  oscillator 
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causes  an  undesirable  frequency  modulation  or  instability  of  the  output 
frequency.  The  instability  range  is  approximately  35  kHz  and,  as  seen 
from  Fig.  B-2,  this  limits  the  achievable  suppression  to  a point  well 
below  the  design  range  of  the  clutter  suppressor. 

It  is  not  surprising  that  the  YIG  oscillator  has  large  frequency 
instabilities  since  it  is  a relatively  inexpensive  source  running  open-loop. 
The  source  can  be  stabilized  by  u.-  ing  it  in  a feedback  loop  and  frequency 
locking  to  a stable  reference.  Such  an  arrangement  is  shown  in  Fig.  4.5. 

The  YIG  signal  is  downconverted  and  filtered  before  being  mixed  with  the 
signal  from  the  reference  synthesizer.  The  signal  out  of  this  mixer  is 
at  10.7  MHz  if  the  YIG  and  the  reference  are  at  the  same  frequency.  The 
frequency  departure  from  10.7  MHz  is  equal  to  the  frequency  error  of 
the  YIG.  This  frequency  error  is  measured  by  the  crystal  discriminator, 
and  the  output  is  used  to  drive  logic  which  applies  a correction  signal 
to  the  sweep  generator.  The  sweep  generator  then  drives  the  YIG  so  as 
to  correct  the  frequency  error. 

Most  of  the  elements  of  this  frequency  lock  loop  are  well  understood. 
Only  the  stable  reference  synthesizer  requires  explanation.  Commercial 
frequency  synthesizers  that  use  direct  digital  synthesis  are  available. 
However,  there  is  none  that  could  be  programmed  to  run  at  the  rates 
required  here.  So  we  designed  such  a unit.  This  is  probably  the  first 
direct  digital  synthesizer  specifically  designed  as  a linear  frequency 
sweeper.  A block  diagram  of  the  design  is  shown  in  Fig.  4.6  where  most 
of  the  lines  are  actually  digital  buses.  The  numbers  by  the  lines 
indicate  the  number  of  bits  in  that  bus. 

At  the  start  of  the  sweep,  the  value  in  the  start  register  is  on 
the  tri-state  bus  and  goes  through  the  summer  to  the  frequency  register. 

Then  the  contents  of  the  sweep  register  are  on  the  tri-state  bus  and 
remain  there  until  the  end  of  the  sweep.  This  bus  value  is  then  added 
to  the  contents  of  the  frequency  register,  and  the  results  are  stored  in  the 
frequency  registei.  This  continues  until  the  contents  of  the  frequency 
register  exceeds  the  value  in  the  end  register.  Then  everything  is 
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reset  and  the  sweep  starts  over.  Clearly,  the  contents  of  the  frequency 
register  increases  linearly  in  time.  The  linear  frequency  is  changed  to 
quadratic  phase  by  a similar  process,  i.e.,  contents  of  the  frequency 
register  are  added  to  the  contents  of  the  phase  register  and  the  result 
is  stored  in  the  phase  register.  Since  a sinusoid  is  a modulo  function 
of  phase,  the  phase  register  will  overflow  in  normal  operation.  This 
phase  value  is  converted  to  a sine  value  by  a read-only  memory  programmed 
with  one-half  cycle  of  a sine  wave.  This  digital  sine  value  is  then 
converted  to  an  analog  signal  which  is  low-pass  filtered  to  obtain  the 
swept  frequency  output. 

The  output  frequency  range  is  0 to  2.1  MHz.  At  a sweep  rate  of 

0.67  MHz/ms,  the  frequency  step  size  is  80  Hz.  The  basic  clock  rate  is 
23 

2 Hz  or  8.388608  MHz  and  the  main  advantage  of  direct  digital  synthesis 
is  that  the  clock  stability  alone  determines  the  stability  of  the  synthesized 
waveform. 

Three  things  will  limit  the  stability  of  the  transmitted  signal. 

The  first  thing  is  the  stability  of  the  stable  reference  synthesizer 
which  depends  only  on  the  clock  rate.  There  will  be  no  difficulty 
eliminating  this  problem.  Second,  the  stability  of  the  stable  oscillator 
will  add  directly  to  the  frequency  jitter.  However,  a commercial  unit 
has  been  identified  which  has  the  required  stability.  Third,  the  YIC 
oscillator  must  be  tightly  locked  to  the  stable  reference.  Because  of  a 
lack  of  an  analytic  model  for  the  YIG  oscillator,  the  behavior  of  the 
frequency  locked  loop  cannot  be  predicted  and  must  be  studied  experimentally. 

Any  solid  conclusions  on  low  angle  operation  must  wait  until  the 
clutter  suppressor  and  frequency  stabilizer  have  been  completed  and  tested 
in  the  radar  system.  However,  we  hope  to  cancel  ground  clutter  by 
30  dB  and  operate  at  elevation  angles  less  than  10°. 

5.  Other  Activities 

During  the  course  of  the  experiment,  certain  other  activities  were 
carried  forward  on  a small  scale  at  a low  priority.  Two  of  these  activities 
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were  in  areas  that  have  direct  application  for  low  level  remote  sensing 
at  airports:  low-level  temperature  profile  measurement  and  wake  vortex 
detection.  The  third  activity,  remote  operation  over  phone  lines,  is 
essential  for  any  operational  application.  This  section  discusses  these 
three  activities. 

5.1  Radio  Acoustic  Sounding  System  (RASS) 

Acoustic  signals  propagate  at  a speed  proportional  to  the  square  root 
of  the  absolute  temperature  of  the  air.  So,  if  a Doppler  radar  can  detect 
an  electromagnetic  signal  reflected  from  an  acoustic  signal,  the  temperature 
of  the  medium  can  be  determined  from  the  measured  Doppler  frequency. 

This  technique  is  called  RASS  and  has  been  verified  experimentally  for 
radar  wavelengths  of  8 m (Marshall  et  al. . 1972).  Since  this  wavelength 
is  80  times  the  wavelength  of  the  FM-CW  radar,  it  was  not  known  if  the 
RASS  technique  would  work  using  a 6 KHz  acoustic  source  and  the  10  cm 
radar. 

Previous  work  using  microwave  frequencies  had  failed  because  of  the 
large  acoustic  absorption  at  6 KHz  acoustic  frequency  (approximately  10 
dB/100  m) . Some  unknowns  were  1)  the  effect  of  sweeping  the  radar  frequency; 
2)  the  effects  of  temperature  gradients  and  turbulence  for  these  short 
wavelengths;  3)  the  effects  of  horizontal  wind.  All  of  these  effects 
would  reduce  the  sensitivity  of  the  system,  but  it  was  hoped  that  the 
FM-CW  radar  would  have  sufficient  sensitivity  to  make  up  for  these  losses. 

We  conducted  some  field  tests  with  an  acoustic  source  adjacent  to  the 
radar  antennas  and  both  the  acoustic  beam  and  the  radar  beam  directed 
vertically.  We  were  able  to  track  the  acoustic  pulse  to  about  a 400  m 
height.  There  were  large  and  fast  fluctuations  in  the  return  signal 
strength.  This  was  probably  due  to  small  scale  turbulence  destroying 
the  focus  of  the  acoustic  signal.  Another  effect  was  that  the  return 
seemed  to  be  very  height  dependent.  This  is  probably  because  the  horizontal 


wind  advects  the  acoustic  signal  away  from  the  radar  beam.  Or  it  could 
also  be  that  temperature  gradients  are  changing  the  acoustic  wavelength. 

The  RASS  experiments  lasted  for  three  days  and  were  designed  only 
to  demonstrate  feasibility  of  the  technique  at  10  cm  radio  wavelengths. 

It  was  shown  that  it  is  not  difficult  to  operate  in  a RASS  mode  to 
heights  of  400  m.  With  some  optimization,  heights  of  600  m should  be 
possible.  We  have  since  analyzed  the  effects  of  radar  frequency  sweeping 
and  temperature  gradients.  The  main  conclusions  are  that  frequency 
sweeping  causes  only  a small  loss  in  sensitivity,  but  temperature  gradients 
can  cause  a significant  sensitivity  loss.  So  a practical  RASS  system 
must  be  a method  for  matching  radar  and  acoustic  wavelengths  at  different 
heights.  The  main  problem  with  RASS  seems  to  be  that  horizontal  winds 
blow  the  acoustic  signal  out  of  the  radar  beam.  This  problem  could  be 
solved  by  moving  the  acoustic  source  upwind  and  letting  the  acoustic 
signal  advect  into  the  radar  beam. 

The  RASS  work  consisted  of  the  three  days  of  experiment  and  some 
ar'.'.ysis  of  the  effects  of  wavelength  changes.  We  demonstrated  the 
t -.-Ability  of  low-level  temperature  profile  measurements  when  the 
horizontal  wind  is  not  too  great.  This  may  be  useful  for  predicting  the 
onset  and  dispersal  of  fog  at  airports. 

5.2  Wake  Vortex  Detection 

The  Naval  Oceanic  Systems  Command  (NOSC)  in  San  Diego  has  an  FM-CW 
radar  somewhat  like  ours  but  without  Doppler  processing  capability.  Some 
unpublished  work  reports  that  they  placed  their  radar  near  a commercial 
airport  and  observed  some  very  strong,  unexplained  echoes.  They  decided 
that  one  possible  explanation  was  the  echoes  were  associated  with  the 
wake  vortices  of  the  commercial  jetliners.  They  repeated  the  experiments 
using  a smaller  Navy  plane  that  could  fly  either  by  jet  or  prop,  but 
could  not  detect  any  echoes  associated  with  the  wake  vortex. 
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While  our  FM-CW  radar  was  at  Buckley  Field,  we  arranged  for  some  of 
the  jet  trainers  to  fly  over  the  antennas  at  about  400  m height  and  at 
a speed  of  about  360  mph.  We  could  see  no  evidence  of  returns  associated 
with  the  wake  vortex  of  any  of  the  aircraft. 

The  positive  results  at  San  Diego  were  obtained  with  large  aircraft 
either  landing  or  taking  off.  The  negative  results  we  obtained  and  the 
negative  results  obtained  at  San  Diego  were  with  smaller  military  aircraft 
(prop  and  jet)  in  level  flight.  The  explanation  is  probably  that  only 
the  larger  aircraft  generate  a wake  vortex  that  lasts  long  enough  to  be 
detected  by  radar. 

5.3  Radar  Operation  Over  Phone  Lines 

The  FM-CW  radar  is  designed  so  that  data  recording  and  antenna 
pointing  are  controlled  by  computer.  The  computer  is  located  in  the  radar 
trailer  and  access  to  the  computer  is  through  a standard  keyboard  in 
the  trailer.  The  keyboard  is  connected  to  the  computer  by  an  RS-232 
interface.  So  with  a commercial  modem  for  an  RS-232  interface  and  a 
teletype  designed  for  phone  line  operation,  we  should  be  able  to  control  the 
radar  over  the  phone  lines. 

We  were  unable  to  have  a phone  line  installed  until  the  radar  was 
at  Buckley  Field.  But  then,  with  only  a small  amount  of  effort,  we  were 
able  to  operate  over  phone  lines.  So,  by  remote  control  we  can  measure 
transmitted  power,  point  antennas,  start  and  stop  recording,  monitor  data 
being  recorded,  and  measure  wind  and  other  moments  of  the  spectra. 

6.  Conclusions 

The  conclusions  we  draw  from  the  experience  gained  in  this  experiment 

2 

relate  to  the  radar  operation  and  to  the  C measurements.  First,  are  con- 

n 

elusions  relating  to  the  radar  and  its  capabilities. 
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(1)  Our  basic  design  for  the  FM-CW  Doppler  radar  is  such  that  off-the-shelf 
components  can  be  used  to  assemble  a very  reliable  wind  sensing  radar. 
We  have  had  problems  with  the  power  transmitter  and  the  computer, 

and  these  problems  have  both  been  traced  to  excessive  heat.  The 
power  transmitter  is  similar  to  an  uplink  satellite  transmitter  and 
the  traveling-wave- tubes  used  in  those  applications  can  achieve 
five  years  of  continuous  operation.  We  believe  that  our  radar  can 
be  designed  to  meet  reliability  requirements  of  an  operational 
radar. 

(2)  We  can  measure  winds  to  heights  of  2 km  almost  all  of  the  time  in  the 
spring,  summer,  and  early  fall.  In  the  late  fall  and  winter  we  can 
measure  winds  to  heights  of  500  m almost  all  of  the  time.  Further, 
the  radar  is  an  all-weather  device  that  can  operate  in  clear-air, 
fog,  rain,  or  snow. 

(3)  The  ability  to  operate  at  ranges  of  5 km  along  the  glide  path  depends 
more  on  the  clutter  suppression  capability  of  the  radar  than  on  its 
sensitivity.  The  radar,  as  now  configured,  is  sensitive  enough 

to  see  to  ranges  in  excess  of  5 km  at  glide-slope  angles  if  the 
ground  clutter  does  not  saturate  the  signal  processor. 

2 

Second  are  conclusions  relating  to  . 

2 

(4)  C can  be  considered  as  a single  parameter  lognormal  random  variable. 

n 222  -2/3 

Only  the  mean  of  log  C (where  log  C =0  implies  C = 1 m ) 

n n 2n 

need  be  specified.  The  standard  deviation  of  log  is  about  0.6. 

2 

(5)  As  expected,  the  mean  C values  are  highest  in  the  summer  and 

n 

lowest  in  the  winter. 
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(6)  The  atmospheric  diurnal  variations  are  such  that  C is  at  a maximum 

n 2 

in  the  early  afternoon  and  at  midnight.  The  minimum  C values  are 

n 

generally  seen  at  sunrise. 

2 

(7)  Because  of  the  low  humidities  in  Colorado,  the  C values  for  the 

n 

central  U.S.  and  east  coast  will  be  greater  than  those  reported  here. 
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APPENDIX  A. 


A SIMPLIFIED  EXPLANATION  FOR  PROCESSING  OF  FM-CW 
DOPPLER  RADAR  SIGNALS  FROM  DISTRIBUTED  TARGETS 

R.  B.  Chadwick  and  R.  G.  Strauch 
NOAA/ERL/Wave  Propagation  Laboratory 
Boulder,  Colorado  80302 

ABSTRACT 

The  concept  of  FM-CW  Doppler  processing  is  explained  by 
showing  it  to  be  a special  case  of  a correlation  receiver  with 
a frequency  offset  reference.  For  a linear  FM  signal,  time  delay 
and  frequency  shift  can  be  arranged  to  cancel  each  other,  and  this 
results  in  considerable  simplification  of  the  correlation  receiver. 
The  number  of  range  cells  and  the  number  of  frequency  bins  per  range 
cell  are  shown  to  be  related  to  the  sampling  obtainable  by  a digital 
spectrum  analyzer. 


This  appendix  is  to  be  published  in  IEEE  Trans,  on  Aerospace  and  Electronic 
Systems,  .January  1979. 
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A SIMPLIFIED  EXPLANATION  FOR  PROCESSING  OF  FM-CW 
DOPPLER  RADAR  SIGNALS  FROM  DISTRIBUTED  TARGETS 

R.  B.  Chadwick  and  R.  G.  Strauch 
NOAA/ERL/Wave  Propagation  Laboratory 
Boulder,  Colorado  80302 

I.  INTRODUCTION 

The  simplicity  and  flexibility  of  FM-CW  Doppler  radar  processing 
makes  it  an  attractive  technique.  The  use  of  this  technique  for  over- 
the-horizon  HF  radars  is  discussed  by  Barrick  (1973).  Strauch  et  al . 
(1976)  used  the  technique  at  S-band  with  an  off-line  digital  Fast  Fourier 
Transform  (FFT)  to  measure  Doppler  shifts  associated  with  returns  from 
refractive- index  inhomogeneities  in  the  clear  atmosphere.  Chadwick  et 
al.  (1976)  used  a standard  spectrum  analyzer  to  measure  these  same 
Doppler  shifts  and  showed  the  feasibility  of  real  time  wind  measure- 
ments. The  technique  has  always  been  difficult  to  explain  because  range 
and  velocity  information  are  extracted  and  separated  simultaneously  by 
a single  spectrum  analyzer  or  FnT.  The  purpose  here  is  to  give  a simple 
explanation  by  showing  that  FM-CW  Doppler  processing  is  a special  case 
of  a correlation  receiver.  Also,  the  importance  of  the  number  of  samples 
per  sweep  and  the  number  of  sweeps  coherently  processed  will  be  pointed 
out . 

A different  way  of  explaining  FM-CW  Doppler  processing  is  given  by 
Strauch  (1976). 
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II.  CORRELATION  RllCIilVHR 


Suppose  that  a transmitted  signal,  s(t)  has  an  autocorrelation 

function  R^.(t)  as  shown  in  Fig.  1,  and  that  this  signal  is  used  in  a 

conventional  correlator  as  in  Fig.  2a.  The  range  is  determined  by  the 

reference  delay  t,  and  the  range  resolution,  A,  is  determined  by  the 

width  of  the  autocorrelation  function  envelope.  The  number  of  cycles 

inside  the  envelope  of  the  autocorrelation  function  is  2A/A,  which  is 

2 4 

generally  in  the  range  of  10  to  10  . When  a target  moves  through  the 
range  cell  at  delay  t,  the  output  of  the  integrator  will  trace  out  the 
autocorrelation  function  and  the  frequency  of  this  signal  will  be  the 
Doppler  frequency  associated  with  the  target.  A special  case  of  this  is 
the  widely  used  police  speed  measuring  radar  where  the  transmitted 
signal  and  the  autocorrelation  function  are  sinusoids. 


When  measuring  the  velocity  of  moving  targets  with  a correlation 
receiver,  the  necessary  integration  should  be  accomplished  not  with  a 
low-pass  filter,  but  with  a bandpass  filter,  or  equivalently  a spectrum 
analyzer.  The  configuration  shown  in  Fig.  2b  will  allow  the  measurement 
of  velocity  spectra  associated  with  either  point  or  distributed  targets. 
However,  velocity  sense  (incoming  or  outgoing)  cannot  be  identified  since 
only  real  signals  are  used.  However,  velocity  sense  can  be  obtained  by 
introducing  a frequency  offset,  fQ,  into  the  reference;  this  is  acceptable 
because  the  offset  simply  translates  the  spectrum.  It  makes  no  difference 
whether  the  frequency  offset  is  on  the  receiver  side  or  the  reference 


side  of  the  correlator  for  both)  A correlation  receiver  with  a frequency 
offset  reference  is  shown  in  Fig.  2c.  Zero  velocity  targets  will  now 
appear  .it  f and  velocity  sense  can  be  determined  by  whether  the  return 
is  above  or  below  f . This  type  of  radar  can  measure  range  and  velocity 
spectra  of  a distributed  target  regardless  of  the  form  of  the  transmitted 
signal.  The  performance  of  this  type  of  radar  when  s(t)  is  a completely 
random  signal  has  been  studied  extensively  and  thoroughly  verified 
experimentally  by  Cooper  and  McGillem  (1967),  Cooper,  McCillem,  and 
Waltman  (1969),  and  applied  to  distributed  targets  by  Chadwick  and 
Cooper  (1972). 

A disadvantage  of  this  type  of  radar  is  that  the  correlator  must  be 
timeshared  (or  parallel  correlators  must  be  used)  to  examine  multiple 
range  cells.  A configuration  using  parallel  correlators  (one  for  each 
range  cell)  is  shown  in  Fig.  3.  The  output  of  the  spectrum  analyzer  of 
the  i^  correlator  is  the  velocity  spectrum  for  the  (point  or  distributed) 
target  at  a range  determined  by  delay  x.  and  this  spectrum  is  centered 
at  f . Up  to  this  point,  the  form  of  s(t)  has  not  been  specified.  If 
s(t)  is  a linear  FM  signal  and  if  the  values  of  f.  and  x.  are  chosen 
correctly,  Fig.  3 can  be  greatly  simplified  as  shown  below. 

III.  LINEAR  FM  SIGNALS  AND  DOPPLER  PROCESSING 

Three  requirements  on  the  parameters  in  Fig.  3 can  greatly  simplify 
the  receiver  structure.  First,  let  s(c)  have  constant  amplitude  and  a 
linear  frequency  modulation  as  in  Fig.  4.  The  solid  line  is  the  instantaneous 
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frequency  of  the  transmitted  signal  and  the  dashed  line  is  the  instantaneous 
frequency  of  the  received  signal  reflected  from  a point  target  at  a 
range  ct/2,  where  c is  propagation  velocity  and  x is  the  round-trip  delay 
time.  Note  that  with  a linear  I’M  signal  the  delay  function  in  the  correla- 
tion receiver  can  be  accomplished  by  a time  delay  or  a frequency  shift. 

The  second  requirement  is  that  the  analyzers  are  gated  off  during 
the  time  intervals  (jT,  jT+x),  j=0,l,2  ...  . This  effectively  reduces  the 
duty  cycle  so  that  the  percentage  loss  of  sensitivity  is  100  x/T.  For 
certain  meteorological  targets,  where  delays  are  on  the  order  of  tens  of 
microseconds  and  sweep  lengths  are  on  the  order  of  milliseconds,  this 
sensitivity  loss  is  negligible. 

The  third  requirement  is  that  the  frequency  offset,  f , for  the  i**1 

branch  of  the  receiver  he  equal  to  i/T.  Then  the  solid  and  dashed  lines 

in  Fig.  4 will  be  col inear,  except  at  the  start  and  end  of  the  sweep 

c i 

for  all  received  signals  whose  delay  is  i/B  or  whose  range  is  So, 

during  the  time  the  analyzers  are  gated  on,  the  signal  at  the  output  of 
the  i**'  frequency  offset  is  just  s(t).  The  frequency  offset  is  equal 
to  the  frequency  shift  corresponding  to  range  or  delay.  This  causes  the 
first  of  two  simplifications.  Since,  during  the  processing  time,  the 
frequency  shifts  due  to  the  delay  elements  are  the  same  as  the  frequency 
offsets,  these  frequency  translations  cancel  so  that  all  of  the  delay 
elements  and  frequency  shift  elements  (in  Fig.  3)  can  be  removed.  The 
second  simplification  then  comes  about  because  each  of  the  parallel 
correlators  has  the  same  reference  and  received  signals  and  so  only  one 
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correlator  is  necessary  to  simultaneously  measure  the  velocity  spectrum 
for  each  range  cell.  The  parallel  correlator  structure  of  Fig.  3 there- 
fore reduces  to  the  structure  in  Fig.  5a.  Because  the  frequency  offsets 
are  different  for  each  delay,  the  spectrum  of  the  correlator  output  is 
translated  differently  for  each  delay  or  range.  Therefore,  the  output 
of  the  single  spectrum  analyzer  has  the  spectra  for  each  range  cell 
displaced  in  frequency  as  in  Fig.  5b. 

The  output  can  be  visualized  further  by  assuming  that  the  return 
is  from  uniformly  distributed,  stationary  ground  clutter.  Then  the  signal 
into  the  spectrum  analyzer  will  be  repeat  exactly  from  one  sweep  to  the 
next.  This  signal  will  then  produce  a discrete  spectrum  with  components 
only  at  multiples  of  1/T  which  are  equal  to  f . . The  location  of  these 
discrete  components  in  Fig.  5(b)  are  zero  velocity  points  for  the  i*^ 
range  cell  and  are  labelled  ZVP^.  There  is  a requirement  on  the  signal 
source  to  insure  that  a stationary  return  produces  a discrete  spectrum. 

The  source  must  meet  a frequency  stability  requirement.  However,  there 
is  no  requirement  on  the  sweep- to- sweep  phase  of  the  signal.  The  sweeps 
do  not  have  to  be  coherently  generated. 

Two  additional  points  should  be  made:  first,  the  velocity  sense 
in  the  closest  bin  can  not  be  determined.  The  reason  for  this  is  that 
the  reference  delay  for  the  closest  bin  is  zero  and  the  frequency  offset 
to  cancel  the  delay  is  zero  and  without  a frequency  offset  reference, 
sense  can  not  be  determined.  If  an  FFT  is  used,  the  velocity  sense  in 
the  last  bin  is  also  lost.  This  problem  can  however,  be  overcome  by 
using  a spectrum  analyzer  that  oversamples  and  has  an  anti-aliasing 
filter. 
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Second,  when  the  Doppler  frequency  becomes  high  enough,  the  return 
appears  in  the  next  higher  or  lower  bin,  i.e.,  it  is  ambiguous.  However, 
even  though  this  occurs,  the  frequency  translation  distance  from  a 
velocity  spectrum  to  the  zero  velocity  point  with  which  it  is  associated 
is  still  the  mean  velocity  at  that  range.  The  ambiguity  function  for 
this  type  of  radar  signal  is  discussed  by  Strauch  and  Cbadwick  (1976). 

The  maximum  unambiguous  velocity  is  derived  below. 

IV.  RANGF,  VF-LOCITY,  AND  DISPLAY  PARAMFTF.RS 

Before  deriving  expressions  for  maximum  range,  maximum  unambiguous 
velocity,  number  of  range  cells  per  display,  and  number  of  Doppler 
frequency  bins  per  range  cell,  it  is  necessary  to  make  some  assumptions 
about  the  spectrum  analyzer.  Although  a fully  analog  analyzer  could  be 
used,  a fully  digital  technique  such  as  an  PFT  or  a hybrid  recirculating 
memory  time  compression  spectrum  analyzer  allows  synchronization  of  the 
sampling  to  the  sweep,  and  this  is  desirable.  Assume  that  the  maximum 
frequency  of  the  analyzer  is  f and  6 is  the  frequency  resolution  (or 
the  frequency  per  bin  of  the  output  display).  So  f^/6  is  the  number  of 
frequency  bins  per  display.  (It  will  be  shown  below  that  6 is  the  Doppler 
frequency  resolution  of  the  system.)  Further,  assume  that  the  input  to 
the  spectrum  analyzer  is  sampled  N times  for  each  sweep  of  the  transmitted 
signal  and  that  there  are  m sweeps  coherently  processed.  Thus,  the 
number  of  time  samples  required  for  a processing  interval  is  Nm  and 
these  samples  result  in  f /6  frequency  bins  in  the  output.  The  ratio  of 
these  two  is  characteristic  of  the  spectrum  analyzer.  If 

no.  of  time  samples  _ Nm6 
^ - no.  of  frequency  bins  " f ’ 
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4 

an  FFT  analyzer  with  Nyquist  sampling  is  characterized  by  p 2.  The 

particular  commercial  hybrid  analyzer  we  use  oversamples , requiring  1500  1 

time  samples  for  500  frequency  bins  and  thus  for  that  analyzer,  p = 3. 

From  Fig.  4,  it  is  seen  that  the  maximum  range,  r^,  is  given  by 


r 

m 


— f = 
2B  m 


ATf 


where  A = c/2B  is  the  range  resolution.  Now  T can  be  replaced  by  the 
number  of  samples  per  sweep  multiplied  by  the  sampling  interval  1/pf  . 
Solving  for  the  number  of  range  cells  per  display,  one  obtains 


r 

jn  _ N range  cells 
A p display 

Also,  it  is  seen  that  dividing  the  number  of  freouency  bins  per  display  by 
the  number  of  range  cells  per  display  gives 

no.  of  frequency  bins  _ 
range  cell 

Therefore,  these  two  important  parameters,  number  of  range  cells  per  display 
and  the  number  of  frequency  bins  per  range  cell,  are  very  simply  related 
to  the  sampling  of  the  spectrum  analyzer  and  can  be  easily  changed.  A 
special  case  is  when  the  spectrum  analyzer  takes  all  of  its  samples  on 
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one  sweep.  For  our.  spect rum  analyzer,  this  means  M 1500,  N'/p  = 500, 
and  m = 1.  So  there  are  500  range  cells  and  only  one  frequency  bin  per 
range  cell;  this  is  the  range-only  processing  case.  The  other  extreme 
of  one  range  cell  with  500  frequency  bins  is  obtained  by  taking  three 
samples  per  sweep  and  coherently  processing  over  500  sweep...  All 
values  between  these  two  extremes  can  be  easily  obtained. 

The  maximum  unambiguous  velocity  is  now  easy  to  find  since  it  is  the 
displacement  from  a zero  velocity  point  by  an  amount  equal  to  half  the 
number  of  frequency  bins  per  range  cell.  The  velocity  of  a target  is 
A f(j/ 2 where  f ^ is  the  Doppler  frequency  equal  to  the  number  of  bins 
displacement  from  the  zero  velocity  point  times  frequency  per  bin.  The 
maximum  unambiguous  di splaccment  is  m/2,  so  that  the  maximum  unambiguous 
velocity  v is 


v 

m 


Tli is  can  be  shown  equal  to  A/4T  as  it  must  be. 

The  velocity  resolution  of  the  system  can  now  be  found  by  dividing 
the  maximum  unambiguous  velocity  by  the  number  of  velocity  cells.  This 
shows  that  the  velocity  resolution  is  6 A/ 2 m/s  or  equivalently,  the 
Doppler  frequency  resolution  is  6 llz. 
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V. 


CONCLUSIONS 


L'M-CW  Doppler  processing  has  two  main  advantages,  simplicity  and 
flexibility.  The  processing  technique  can  be  considered  a special  case 
of  parallel  correlators  with  frequency  offset  references.  The  requirement 
of  linear  I'M  modulation  on  the  transmitted  signal  simplifies  the  parallel 
correlators  to  a single  mixer  and  a single  spectrum  analyzer.  The  maximum 
range  and  range  resolution  are  determined  by  the  frequency  sweep  which 
is  easily  changed  by  a low  level  control  voltage.  The  number  of  range 
cells,  number  of  frequency  bins  per  range  cell,  and  the  maximum  unambiguous 
velocity  is  controlled  by  the  sampling  of  the  spectrum  analyzer  which  is 
also  easy  to  change. 


Alt) 


VI . REFERENCES 


Barrick,  D.  E.,  FM-CW  Radar  Signals  and  Digital  Processing,  NOAA  Technical 
Report  ERL  283-WPL26,  July  1973. 

Strauch,  R.  G.,  W.  C.  Campbell,  R.  B.  Chadwick,  and  K.  P.  Moran,  Micro- 
wave  FM-CW  Doppler  Radar  for  Boundary  Layer  Probing,  Geophysical 
Research  Letters,  Vol . 3,  no.  3,  p 193-196,  March  1976. 

Chadwick,  R.  B.,  K.  P.  Moran,  R.  G.  Strauch,  G.  E.  Morrison,  and  W.  C. 
Campbell,  A New  Radar  for  Measuring  Winds,  Bulletin  of  the  American 
Meteorological  Society,  Vol.  57,  no.  9,  p 1120-1125,  Sept.  1976. 

Cooper,  G.  R.  and  C.  D.  McGillem,  Random  Signal  Radar,  Technical  Report 
EE67-11,  Purdue  University,  Lafayette,  Indiana,  June  1967. 

Cooper,  G.  R.,  C.  D.  McGillem  and  W.  B.  Waltman,  Use  of  Wideband  Stochastic 
Signals  for  Measuring  Range  and  Velocity,  EASCON  '69  Record,  p 305-311, 
1969. 

Chadwick,  R.  B.  and  G.  R.  Cooper,  Measurement  of  Distributed  Targets 

with  the  Random  Signal  Radar,  IEEE  Trans,  on  Aerospace  and  Electronic 
Systems,  Vol.  AES-8,  no.  6,  p 743-750,  Nov.  1972. 

Strauch,  R.  G.  and  R.  B.  Chadwick,  Measurement  Capabilities  of  FM-CW 
Doppler  Radars,  17th  Conference  on  Radar  Meteorology,  American 
Meteorological  Society,  Oct.  1976. 

Strauch,  R.  G.,  Theory  and  Application  of  the  FM-CW  Doppler  Radar,  Ph.D. 
Thesis,  University  of  Colo.,  1976. 


A 1 1 


Figure  1.  Autocorrelation  function  of  transmitted  signal. 


Signal  

Source  “"T  s(t) 


Non  Moving  Target 
Target 


(c) 


Figure  2(a).  Conventional  correlator,  (b)  Bandpass  correlator. 

(c).  Bandpass  correlator  with  frequency  offset  reference. 


A1 2 


0 


U ' 

ZVP]  Z' 


CPI  ■ ( ■■) . " . 


'radar  oubjut  for  five  range  • 1 r, 
rith  a rrooina  target  in  each  cell 


APPENDIX  B.  l-'.oqui rod  Frequency  ''..In  lit;  to  Suppress  Ground  Glutt or 


In  Appendix  A,  it  was  pointed  out  that  the  FM-CW  radar  had  no 
requirement  on  the  phase  of  the  transmitted  signal,  but  there  was  a re- 
quirement on  the  starting  frequency  of  the  FM  sweep.  In  this  section, 
that  requirement  is  quantified,  and  this  allows  us  to  determine  maximum 
clutter  suppression  as  function  of  range  and  instability. 

A simplified  block  diagram  of  the  FM-CW  radar  is  shown  in  Fig.  B-l. 
In  this  analysis,  the  amplitudes  of  the  transmitted  and  received  signal 
are  not  important  and  will  be  taken  as  unity.  The  transmitted  signal 
for  one  sweep  of  length  T is 

0 < t < T 


where  the  bandwidth  is  B Hz  and  the  starting  frequency  is  f Hz.  Assume 
the  return  signal  is  from  a point  ground  clutter  target  and  has  been 
delayed  by  t so  that 


s(t)  = sin  2n<J>T(t) 


V(>  ■ i?  *2  * V 


s(t-t)  « sin2u$T(t-i:) 


111 


and  the  signal  out  of  the  mixer  is  then 


r (t)  = cos2TnJ>r(t) 

♦r(t)  = 4>T(tJ  - <J>T(t-x). 


The  phase  term  at  the  output  of  the  mixer  is  then  shown  to  be 


- j tt 


IB  2 
2 T 


f T 

O 


So,  during  a sweep,  r(t)  is  a sinusoid  with  frequency  pro]>ort iona 1 to 
the  delay  of  the  return  signal,  and  with  a phase  equal  to  delay  times  RF 
sweep  starting  frequency.  Tt  is  clear  that  any  instability  or  jitter  in 
this  starting  frequency  will  cause  r(t)  to  change  from  one  sweep  to  the 
next,  limiting  the  effectiveness  of  any  device  designed  to  cancel  ground 
clutter. 


To  determine  how  much  this  frequency  jitter  will  limit  the  maximum 
clutter  suppression,  we  must  find  the  spectrum  of  r(tj.  For  a single 
point  clutter  target  at  a certain  range,  r(t)  will  be  a sinusoid  and  if 
f has  some  jitter,  this  sinusoid  will  change  phase  randomly  at  time 
T,  2T,  ...  This  is  very  much  like  the  signal  for  a phase  shift  keyed 
(PSK)  digital  communication  system.  Lundquist  (1909)  has  derived  a very 
general  expression  for  the  power  spectrum  of  a PSK  signal  x(t)  where 
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x(t)  = V COS(U)  t + <Mt)) 
c c 

OO 

<Ht)  = l an  g(t-nT)  . 

n-oo 

Here,  we  are  interested  in  the  special  case  where  g(t)  = 1 for  0 < t < T 
and  is  zero  elsewhere.  Also  we  desire  r phase  levels,  equidistant  and 
equ iprobabl e , i . e . , 


a 


(k)  = 


a + 
o 


k-1 

r 


a 


Prob[a(k)]  = I . 


It  will  be  necessary  to  let  r become  infinite  for  this  model  to  approximate 
our  situation.  The  parameter  a will  then  be  the  radian  phase  jitter 
in  r(t)  due  to  the  frequency  jitter  in  fQ. 


The  power  spectrum  of  this  PSK  signal  is 


where 


S(u>c+o>)  = f- 


2tt 


sin(  j1-) 

~^r 


P A sin(lj 
sir,(yr-) 


Note  that  there  is  a continuous  part  and  a discrete  part  of  this  spectrum. 
The  continuous  part  is  due  to  the  changes  in  phase  of  r(t),  which  arc  caused 
by  the  jitter  in  sweep  starting  frequency,  f . So,  the  power  in  the 


continuous  spectrum  can  be  called  the  jitter  power,  J,  and  the  power  in 
the  discrete  spectrum  can  be  called  the  signal  power,  S.  The  signal 
power  to  jitter  power  ratio  is  then 


S 

J 


To  evaluate  this  limit,  note  that  as  r becomes  large,  the  sine  term  in 

2 

the  denominator  of  F can  be  replaced  by  its  argument.  Then  the  r terms 
in  S/J  will  cancel  leaving 


S _ A 
J " 1-A 


where 


s i n a/2 

a/2 


The  range  of  phase  change  of  r(t)  is  then  proportional  to  the  range  of 
starting  frequencies  or  the  frequency  jitter 


A<(>  = 2iTTAf  = a . 

r o 

The  signal-to-j i tter  ratio  can  now  be  plotted  as  a function  of  the  frequency 
jitter  with  delay  as  a parameter.  This  plot  is  Fig.  B-2. 

The  main  use  of  Fig.  B-2  is  to  determine  the  maximum  achievcablc 
clutter  suppression.  Since  a clutter  suppressor  can  only  suppress  signals 
with  a discrete  spectrum  the  si gnal -to- j i tter  ratio  is  the  maximum 
clutter  suppression  that  is  achievable  for  that  particular  x and  Af^.  For 


example,  if  it  is  necessary  to  suppress  ground  clutter  from  450  m range 
by  30  dB,  the  frequency  jitter  must  be  less  than  8 kHz.  The  ground 
clutter  suppressor  that  is  presently  under  construction  has  a suppression 
potential  of  48  dB.  To  achieve  this  for  ranges  of  450  m or  less,  a 
frequency  stability  of  about  1 kHz  is  needed. 

Reference 

I.undquist,  L.  (1969),  Digital  PM  Spectra  by  Transform  Techniques,  BS'I.J, 
vo 1 . 48.  no.  2,  p.  397-411. 
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Simple  block  diagram  of  FM-CW  radar. 
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Figure  B-l. 


Signal-to-J itter  Ratio  (dB) 


Sweep  Starting  Frequency  Jitter  (kHz) 

Figure  B-P..  Frequency  stability  required  for  a given  signal- to- jit  ter  ratio. 
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AtTI.NDIX  C.  Comparison  of  Clutter  Power  for  a Pulse  Radar  and  a CW  Radar 


The  wind  shear  problem  requires  operation  at  low  elevation  angles 
and  short  ranges.  This  appendix  compares  how  easily  two  different  radars 
can  operate  at  low  angles  and  short  ranges  by  finding  the  ratio  of  return 
clutter  powers.  The  critical  aspect  of  operation  at  low  elevation  angles 
and  short  ranges  is  whether  or  not  the  receiver  saturates.  During  the 
time  the  receiver  is  saturated  by  ground  clutter,  wind  measurements  are 
not  possible.  If  the  radars  have  similar  first  stages  in  the  receivers, 
the  results  here  tell  which  is  the  most  susceptible  to  saturation. 

Assume  the  geometry  shown  in  Fig.  C-l  where  the  antenna  is  h meters 
above  the  ground  and  looking  at  an  elevation  angle  e.  Let  be  the  angle 
from  antenna  boresite  and  assume  the  antenna  beam  is 

G(<p)  = Cj  |*|  < c 

= Cj/V72  |*|  > e . 

The  shape  of  the  beam  for  angles  less  than  the  elevation  angle  is  not 
important  because  this  portion  of  the  beam  does  not  see  the  ground. 

Assume  that  the  ground  has  radar  cross  section  per  unit  area  of  o° . 
The  incremental  power  returned  from  a incremental  range  element,  dR  at 
range  R,  in  a plan  view  sector  0 radians  wide  is 


Cl 


r * 


pt  n2  (4>j  2 ° 

dP  = — = — : — OR  dR 

c 3 4 

(4ttJ  R 


where  I*  is  the  transmitted  peak  power.  Since  this  is  a comparison  where 
a ratio  of  powers  is  formed,  the  common  terms  will  cancel  out.  So  rather 
than  carry  them  along,  we  will  drop  them  when  they  become  multiplicative 
factors.  The  total  clutter  power  can  then  be  found  by  integrating  over 
the  illuminated  range  interval. 


,,2  , , , 

P = PT  / 2 (- dR 
C 1 R j R3 


where  terms  that  will  eventually  cancel  have  been  dropped. 


Because  of  the  low  elevation  angle. 


♦ ~ 


eR+h 


and  the  clutter  power  expression  becomes 


R, 


P » P / 2 

C T Rj  (eR+h)3 


This  integral  evaluates  to 


P-r 


(eR  j +h) 2 (cR2+h)2 


C2 


lor  a CW  radar  the  maximum  illuminated  range  is  infinite  so  the  second 
term  goes  to  zero.  The  minimum  illuminated  range  can  be  taken  as  the 
smaller  of  the  far-field  boundary  and  the  point  at  which  the  antenna 
beams  intersect.  Mere  the  far-field  boundary  will  be  used.  So  the 
minimum  illuminated  range  for  the  CW  radar  will  be  taken  as 

r|  = 1)2/A  . 

Then,  the  ratio  of  clutter  power  for  a pulse  radar  to  that  for  a CW  radar 
is  given  by, 

jz  pT  (YctyhV  /cvhV 

pc  ’ ptLvRi+7  ’ v2*hJ  . 

where  primed  terms  arc  for  the  CW  radar. 

This  ratio  does  not  change  greatly  when  e and  h are  changed  over 
reasonable  ranges  (but,  of  course,  the  characteristics  of  a single  specific 
radar  do  change  as  c and  h are  changed) . The  clutter  power  of  a CW  radar 
is  independent  of  range  cell  while  for  a pulse  radar  the  clutter  power  is 
approx imately  proportional  to  the  range  cell  size.  So  this  ratio  varies 
with  the  range  cell  size.  A plot  of  pulse  clutter  power  over  CW  clutter 
power  is  shown  in  Fig.  C-2  for  some  reasonable  parameters.  Whenever  this 
ratio  is  positive  in  dB,  the  pulse  clutter  problem  is  more  severe  than 
the  CW  clutter  problem. 

C3 


This  analysis  is  simplified  because  of  the  assumption  of  uniform 
reflectivity  of  the  ground.  Terrain  with  man-made  objects  is  difficult 
to  model  and  usually  experimental  techniques  must  be  used.  However,  a 
conclusion  of  this  appendix  is  that  if  the  ground  close  to  the  radar  is 
relatively  smooth  and  free  of  man-made  objects,  an  l-'M-CW  radar  will  be 
easier  to  operate  at  short  ranges  and  low  elevation  angles. 


Figure  C-l.  Geometry  for  ground,  clutter  comj>arinon. 


o o o o o 
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Figure  C-2.  Ratio  of  clutter  pcwer  (pulse/CW)  for  vari<  w •••  tenna  H ”■  fc<  re. 
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APPENDIX  D.  Histograms  of  C 

n 

This  appendix  gives  probability  histograms  of  the  measured  C ^ 

values  for  each  of  the  nine  height  bins.  The  histograms  are  presented 

two  different  ways;  first  on  a monthly  basis  and  then  on  an  hour  of 

day  basis.  The  numbers  represent  the  percentage  of  time  the  measured 
2 

was  within  that  bin.  A number  of  .000  indicates  the  percentage 
of  time  was  less  than  .0005,  but  greater  than  zero.  A blank  indicates 
that  the  percentage  of  time  was  zero. 

Por  the  last  four  months  of  the  experiment,  the  radar  was  located 
in  an  urban  area  where  the  ground  clutter  was  large.  So,  the  lower  two 
height  intervals  could  not  be  used. 

The  hour-of-day  (MST)  histograms  are  presented  in  a similar  format, 
except  only  one  height  bin,  805  m,  is  included.  The  data  obtained  in 
the  last  four  months  of  the  experiment  were  not  used  for  the  hour-of- 
day  histograms. 

A more  convenient  way  of  looking  at  this  data  is  used  in  Section  3. 
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TABU:  11.8  Monthly  probability  C ~ histograms  at  height  1288 


TABLh  0.9  Monthly  probability  C “ histograms  at  height  1449 
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